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1 Introduction 



The application of perturbative qcdEI is not straightforward even for reactions 
with large momentum transfer Q. In higher order corrections soft and collinear 
momentum regions may lead to large contributions of order (a s \og(Q / m q )) n 
where m q denotes a light quark mass. These terms give 0(1) corrections which 
destroy the validity of the perturbative treatment. The applications of pertur- 
bative QCD are limited to phenomena where such terms arc either cancelled 
or can be controlled with improved treatment (resummation). Fortunately, 
in periurbative QCD the soft and collinear structure is relatively well under- 
stood □. The main features are summarized by fundamental cancellation and 
factorization theorems valid in all orders in perturbation theory. 

This lecture will show explicitly how the cancellation and factorization the- 
orems "work" in next-to-leading order applications. I describe how to combine 
virtual and real next-to-lcading order infrared singular amplitudes into finite 
physical cross sections. The method of calculating amplitudes in leading, and 
next-to-leading order have been described in great detail by Lance DixonB. He 
also discussed the asymptotic properties of amplitudes in the soft and collinear 
regions as well as the soft and collinear singular terms appearing in one-loop 
amplitudes. I describe the singularities as they appear in various cross sec- 
tion contributions (loop contributions and bremsstrahlung-contributions). I 
consider their universal features and show explicitly that they cancel for in- 
frared safe quantities. The examples will always be chosen from the physics 
of jet-production. In order to build a numerical program for efficiently calcu- 
lating hard scattering cross sections in next-to-leading order (NLO) accuracy 
one should use a well defined algorithm for the analytic cancellation of the soft 
and collinear singularities. Several methods are used in the literature. I shall 
focus on the subtraction method. 



2 Cancellation and factorization theorems 

2. 1 KLN cancellation theorem 

In simple inclusive reactions, such as the total cross section of e + e~ annihi- 
lation into quarks and gluons, the soft and collinear contributions cancel □. 
That is a consequence of the KLN theorem. In the simplest examples, only 
one high momentum transfer scale is relevant, the effective coupling becomes 
small and the cross section can reliably be calculated in power series of the 
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effective coupling up to small power corrections 
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as is the running coupling constant, [i is the ^normalization scale, B2 is a 
known constant given by the NNLO calculation!] and 60 is the first coefficient 
in the beta-function 

where nj denotes the number of quark flavours. The truncated series is /1- 
dependent but the fi dependence is 0(as^ n+1 ') if the cross section is calculated 
to 0{a s n ). 

The KLN theorem remains valid also for integrating over final states in 
a limited phase space region, as is the case of jet production. The Sterman- 
Weinberg two-jet cross sectionQ is defined by requiring that all the final state 
partons are within a back-to-back cone of size 5 provided their energy is less 
than e^/s. At NLO 



c 2 jct = cr s -w(s,e,S) 

= otot - 0gsl( al1 E > eVs, all % > S) 
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where cro = 47ra 2 /3s. The cancellation theorem here is hidden in the calcula- 
tion of the total annihilation cross section. 



2.2 Factorization theorem 

The initial state collinear singularities, which in general do not cancel, are uni- 
versal and process-independent in all orders in perturbation theorya. Therefore 
they can be cancelled by universal collinear counter terms generated by the 
'renormalization'of the incoming parton densities. The rule for defining the 
finite part of this counter term is fixed by the factorization scheme. As in the 
case of ultraviolet renormalizationl3, the physics is unchanged under a change 
of the factorization scheme, provided the parton densities are also changed 
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suitably. This feature is expressed by the Altarelli-Parisi evolution equation 
of parton densities. The collinear subtraction terms define the kernels of the 
evolution equations. 

The differential cross section for hadron collisions can be written as 

d(JAB(PA,PB) = J dx 1 dx 2 f a /A(xA)fb/B{xB)dd- ab (xAPA,XBPB) , (5) 
ab 

where A and B are the incoming hadrons, pa and pb their momentum, and 
the sum runs over all the parton flavours which give a non-trivial contribution. 
The quantities da ab are the subtracted partonic cross sections, in which the 
singularities due to collinear emission of massless partons from the incoming 
partons have been cancelled by some suitable counter terms. 

According to the factorization theorem, the subtracted cross section is ob- 
tained by adding the collinear counter terms to the unsubtracted cross section. 
The latter quantity can be directly calculated in perturbative QCD. Due to 
universality, eq. (|^) applies also when the incoming hadrons are formally sub- 
stituted for partons. In this case, we are also able to evaluate the partonic 
densities, which at NLO read 

f a /d(x) = S ad S(l - x) - |£ (^P a/d (x, 0) - K a/d (x)^j + O {as 2 ) , (6) 

where P a / d (x,0) are the Altarelli-Parisi kernels in four dimensions (since we 
will usually work in 4 — 2e dimensions, the in the argument of P a / d stands for 
e = 0) and the functions K a / d depend upon the subtraction scheme in which 
the calculation is carried out. For MS, K a i d = 0. Writing the perturbative 
expansion of the unsubtracted and subtracted partonic cross sections at next- 
to-leading order as 

da ab = da a ° b ] + da£ , da ab = da a ° b ] + da^ , (7) 

where the superscript 0(1) denotes the leading (next-to- leading) order contri- 
bution, we have 

d^ b \ Pl ,p 2 ) = da { ° b \ Pl , P2 ) (8) 
d^ipun) = da a 1 b ) ( PllP2 ) + da c a °r t (p 1 ,p 2 ) (9) 
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where 



da c a °r t (Pi,P2) = ^T,J dx {j P d/ a {x, 0) - K d/a (x)j da^(x Pl ,p 2 ) 

+ ff J dx (J= p d/b(%,0) - K d/b (x)\ da a °](pi,xp 2 ) . 

(10) 

Eq. (|o|) defines the collinear counter terms for any finite hard scattering cross 
section for processes with quarks and/or gluons in the initial state. Notice that 
in this equation the Born terms da^ are evaluated in 4 — 2e dimensions. 

The KLN theorem and the factorization theorems, constitute the theoret- 
ical basis of the description of scattering processes of hadrons in perturbative 
QCD. Those physical quantities for which these theorems remain valid are 
called infrared safe^. These theorems constitute the necessary consistency con- 
dition for the validity of the fundamental assumption of the QCD improved 
parton model. This assumption is that for the case of infrared safe quanti- 
ties the perturbative QCD predictions given in terms of partons are a good 
approximation to the same quantities measured in terms of hadrons (up to 
power corrections which are small at high momentum scales). 

Provided the higher order corrections at a given order are larger than 
the power corrections, one can systematically improve the accuracy of the 
predictions by calculating terms of higher and higher order. Indeed the analysis 
of the experimental results required the inclusion of higher order radiative 
corrections for a large number of measured quantities. 

3 Jet cross sections at next-to-leading order 

We consider the cross sections for three jet production in e + e~ annihilation 
and two-jet production in hadron-hadron collisions. These cross sections are 
proportional to at least one power of as and are studied experimentally with 
high precision. 

3.1 Three-jet production in e + e~ annihilation 
Let us consider the process 

e-(k-)+e + (k+) =o 1 (pi) + ... +a n ( Pn ) (11) 

a In other words a measurable is infrared safe if it is insensitive to collinear splittings of 
partons and/or emission of soft gluons. 
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where on denote quarks or gluons and n = 3,4. The amplitudes of these pro- 
cesses A^ 1 '^ are known in the tree approximation [i = 0) for n — 3, 4, 5 partons 
and in the one-loop approximation (i = 1) for (n — 3) partons. Lance Dixon 
explained how to calculate these amplitudes quickly with modern techniques. 
It is convenient to consider the squared amplitude divided by the flux and the 
spin averaging factor 8s (s = 2fc + fc_) : 

colour 

spin 



and 



colour 

spin 

where {w;} m ,m+n is a short-hand notation for the list of variables 

and we indicated the flavour and momentum dependence only for the ip func- 
tions. The one-loop corrections to the production of three partons are given 
by V (3 ' 1} - They were calculated first by R.K. Ellis, Ross and Terrano (ERT)la. 
Recently, a new derivation using the helicity method, where the orientation 
with respect to the beam direction is not averaged, was given by Giele and 
GloverEI. 

The physical cross sections are obtained by integrating the product of the 
■0 functions and some "measurement functions" Sx over the corresponding 
phase space volume 

dcr nl ° = da Bom + da vilt + dcr rcal , (14) 



where 



da Bo ™(s,X) = Y, d°-^3({«}i,3^)# 3 ({B}u) (15) 

{ai}l,3 

da vM (s,X) - J2 ^e- S xA{Pih^X)dM{PiK 3 ) (16) 
da™\s,X) = Y ^f e -Sx A ({Pi}iA;X)dM{PiK 3 ) (17) 

{ a l}lA 
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where X stands for the measured physical quantity and 



1 f " rid— i n 

*>» = n (2^( 2 -)^ (d) ( fc + + k - x>) 

with d = 4 — 2e. As a result of complete flavour sum the final particles behave 
as though they were identical; this explains the identical particle factor of nl. 
All quantities are calculated in d = 4 — 2e dimensions. The singular terms 
appear as single or double poles in e. The singularities, however, cancel in the 
sum (Q). 

Infrared safe measurement function. The cancellation of the soft and 
collinear singularities of the virtual corrections against the singular part of 
the real contribution is independent of the form of the measurement functions 
provided they are insensitive to collinear splitting and soft emission. This 
means that one obtains the same measured result whether or not a parton 
splits into two collinear partons and whether or not one parton emits another 
parton that carries infinitesimal momentum. A physical quantity that is de- 
signed to look at short distance physics should have this property, otherwise 
it will be sensitive to the details of parton shower development and hadroniza- 
tion. The mathematical requirements for iSx,3 and Sx,4 are that <Sx,4 should 
reduce to <Sx,3 when two of the outgoing partons become collinear: 

SxM,rf,(l- A)p£,Ap£) =S x , 3 (rf,rf,p$) (19) 

for < A < I plus similar conditions where the A and 1 — A factors are inserted 
to all possible pairs of the momenta in Sx,4,(Pi i P% > P3 , P%) ■ 

Example of an infrared safe observable. For the sake of illustration I give 
the definition of the measurement function for the shape variable thrust T 

n 

S T ,n = fi(T-T n (pf,p% ) , where r„=max^ . (20) 

i=l 

Thrust is well defined for arbitrary number of final-state particles. It is easy 
to check that it satisfies the conditions of infrared safety formulated with the 
help of eq. @. 

Thrust measures the sum of the lengths of the longitudinal momenta of 
the final particles relative to the thrust axis n chosen to maximize the sum. 
For two-particle final states its value is 1. Its allowed range changes with the 
particle number, therefore its differential distribution is only well defined after 
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particles 


T = 1 


Discontinuous 


Three 


particles 


1 > T > 2/3 


in particle 


Four 


particles 


1 > T > 1/ a/3 


multiplicity 


oo 


particles 


1 > T > 1/2 





Table 1: Range of thrust for various numbers of partons 



some smearing]]. Carrying out the phase-space integral for the Born term one 
gets 



1 dcr Born a, 4 



Co dT 2ir 3 



2(3T 2 -3T + 2) , (2T-l\ 3(3T-2)(2-T) 
In 1 ' 



(1 — T)T \1-T J (1 -T) 

(21) 

In the classic paper of ERT the singular pieces have been evaluated an- 
alytically and it was demonstrated that indeed for infrared safe quantities 
they cancel as is required by the KLN theorem. The remaining finite next-to- 
leading order cross section is suitable for numerical evaluation. ERT calculated 
the distribution for the shape variable C. But their analytic result can be used 
to calculate any infrared safe three-jet-like quantity in next-to-leading order. 
It is convenient to give the one-dimensional distributions of various three jet 
measures generically denoted as X = t,C...Ej in a form that satisfies the 
renormalization group equation 

1 d ° asifl) -A x (X) + ( \A x (X)2,b lo g ^/S)+B x (X)] (22) 



do dX 2n \ 7r 

A(x) and B(x) are scale-independent functions. Their values are tabulated 
for many quantities in ref. EEI. The next-to-leading order expression is scale 
independent up to 0(as 3 ) 



d f da 
djj? [dX 



= 0(a s 3 ) . (23) 



The size and sign of the corrections is rather different for the various jet mea- 
sures. In many cases the corrections are substantial (rs 30%) even at LEP 
energy. Thanks to the technical development described by Lance Dixon, the 
NLO calculation will soon be available also for four-jet production ( ip^ X \). 



^Smearing is required also by hadronization effects. The typical width of a Gaussian 
smearing is AT = m^/yfs s=s 0.02 . 



8 



3.2 Jet cross section in hadron-hadron collisions 

At the Tcvatron, multijet cross sections are observed up to six jets 0. The 
analysis of the data requires the evaluation of the amplitudes of the parton 
processes 

Oi(pi) + a 2 (p 2 ) -> a 3 (p 3 ) + ... + a n (p n ) (24) 

where aj denotes parton flavour labels, are their four-momenta and n is the 
number of the participating partons. The data are rather precise for two- and 
three-jet like quantities therefore the comparison with the theory has to be 
done at next-to-leading order. For this purpose the tree amplitudes have to 
be known for n = 4, 5, 6 while in the case of n = 4, 5 we have to know also the 
one-loop amplitudes. It is convenient again to introduce ip functions giving the 
squared amplitude divided by the flux and spin averaging factors 

^"• 0) ({a*}i,«;Mi,n) = 9 , \ f ; E \ A ^ 0) ( 25 ) 

2s oj(ai)ui(a2) I 



spin 

l 



2s <jj(a\)u){a2) 

J2 (-4 ( "' 0) A {n ' 1] * + A (nfi) * „4 ( "< 1} ) (26) 



spin 

where s = 2p\p2, A^ n '^ denote the tree- (i = 0) and one- loop (i = 1) am- 
plitudes of the process (|24|) (helicity, flavour and momentum labels are all 
suppressed), and u>(a) is the number of colour and spin degrees of freedom for 
the flavour a, in 4 — 2e dimensions 

tu(q) = 2N c , iu(g) = 2(1 - e)(7V c 2 - 1) . (27) 

The hard-scattering cross section is decomposed into four contributions 

da al ° = da Boin + da vht + dcr real + da count , (28) 

where 

d°%£(pi,P2lX) = E ^° ] Sx^2({Pl}3,n:X)x 
{a t } 3 ,n 

d<f> n -2(pi,P2;{pi}3,n) (29) 

da^ t a2 (p 1 ,p 2 ;X) = E ^ nA) Sx,n~2({pih, n ;X)x 

{a ( } 3 ,n 
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dc/> n -2(Pl,P2;{pi}3,n) (30) 
{o)}l,n+l 

d(/) n ^i(pi,p 2 ;{pi}3,n+i) (31) 

where d4> n -2(j>i,P2] {pi}3,n) is the phase-space volume for n — 2 final parti- 
cles with total energy defined by the incoming four- momenta p\ and pi , X is 
a generic notation of the measured physical quantity, Sx,n denotes the mea- 
surement functions of X defined in terms of n partons and the counter-term 
cross sections dcr^° u (pi,P2i X) are given by eq. ([To|). The loop corrections 
■0(44) have been obtained for the spin- independent case by Ellis and SextonlI3; 
the spirwdependent one-loop corrections have been obtained by the helicity 
method 113 and very recently the one-loop corrections to the five-parton pro- 
cesses ijj^ 5 ' 1 ' have also been calculated (see the lecture of Lance Dixon). All 
four terms contributing to the hard scattering cross section ( p8h have to be cal- 
culated in 4 — 2e dimensions. The individual terms have singular 1/e and 1/e 2 
contributions. The singularities cancel in the sum, provided the measurement 
functions satisfy the conditions of infrared safety formulated in subsection 3. if . 
The hard scattering cross section given by eq. ( |28| ) is finite and, provided the 
cancellation of the singular terms is achieved analytically, it is suitable for the 
numerical evaluation of physical cross sections with the use of eq. (0) . 



4 Methods of analytic cancellation of the singularities 

Although eqs. (14), (p8|) define finite cross-sections, they cannot be used di- 
rectly for numerical evaluation since the singular terms in the real contribu- 
tions are obtained by integrating over the soft and collinear kinematical range. 
Since the phase space is large and its boundary due to the presence of arbi- 
trary measurement functions is complicated, analytic evaluation is impossible. 
Fortunately, in the soft and collinear regions the cross-sections and the mea- 
surement functions have a simple universal behaviour such that the integration 
relevant for the calculation of the singular terms becomes feasible analytically. 
This feature can be implemented in two basically equivalent methods. 

In the first case ( phase space slicing method ) one excludes (slices) from 
the numerical integration domain the singular regions such that the numerical 
integrations becomes well defined. In the excluded regions the t/j functions, 
iS functions and the relevant phase-space factors can be replaced with their 

c In the case of hadron-hadron collisions the measurement functions should also fulfil the 
condition that Sx,n— 1 should reduce to Sx,n— 2 when one of the partons becomes parallel 
to one of the beam momenta. 
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limiting values and the integrals in these regions are performed analytically. 
The boundary of the excluded regions are defined with some small parameters 
(invariant mass parameter or some small angle and energy-fraction parameter 
as in the case of the Sterman- Weinberg jet cross-sction) . Let us consider as 
illustration a one-dimensional problem with integration domain < x < 1 and 
an integrand which has a simple pole at x = 0. One slices the integration 
region into two pieces, < x < 6 and 6 < x < 1. We choose S -€1 1, thus 
allowing us to use the simple approximation F(x) — > -F'(O) for < x < S. This 
gives 



Js % 

Now the second integral can be performed by normal Monte Carlo integration. 
As long as 6 is small, the sum of the r twii±erms will be independent of 8. About 
the first use of this method, see refs.oll3u3. The method was further developed 
using systematically the universal features of the soft and collinear limits by 
Giele, Glover and KosowerBO. The actual implementation is not completely 
straightforward. The numerical integration over the real contribution has to 
be done with a certain accuracy. Since the integrand near the boundaries of 
the singular regions increases steeply the parameters defining the boundaries 
must not be very small. Furthermore, the result of the analytic integration 
is approximate since the integrand is replaced with its limiting value. The 
result is more and more accurate as the parameters defining the boundary of 
the singular regions become smaller and smaller. Fortunately, the measured 
cross sections have finite accuracy, which sets the required precision for the 
theoretical evaluation. In practical applications one compromises over these 
conflicting requirements to achieve the best efficiency. 

The other method is called the subtraction method. This method takes 
account of the fact that the singular behaviour of the real contribution is 
just some simple pole with well known simple residue. After subtraction the 
numerical integration over the subtracted integrand becomes convergent. The 
quantity subtracted should be added back on. The subtraction terms, however, 
are simple, the dimensionally regulated singular integrations can be carried 
out analytically. This can be illustrated again with a simple one-dimensional 
integral. We write 



I 




(32) 



I 
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= [F(x)-F(0)] + -F(0). (33) 

Jo x e 

The integral can now be performed by Monte Carlo integration. 

This method was first used for QCD calculations by R.K. Ellis, Ross, and 
TerrancH and it was further developed, with systematic use of the simplicity 
of the. .soft ..and collinear limit in refEJ. More applications can be found in 
refs. li3'ESE3£y. The disadvantage of this method is that outside the singular 
region the values of the physical parameters at a given phase-space point will 
change if we take the limit corresponding to the soft or collinear configurations 
of the subtraction terms. In the case of Monte Carlo numerical evaluation 
this can lead to relatively large fluctuations in the binned distributions. This 
problem can be avoided by using Gaussian smearing over the bin sizeLS. 



5 Virtual contributions 



In the following I shall consider only the squared matrix elements of pro- 
cess . The virtual contributions to the cross sections are given be eq. (|o|) . 
We are interested in the form of the singular terms of the functions ip^ n,1 \ It 
turns out that they have a very simple general structure: 

\-~ E - - E ^ ( "' 0) ({aih, n , {Plh, n ) 

[ 6 1=1 6 1=1 J 

+ ^ 1) ({ai}i, n ,{Pi}i,n) (34) 
where we introduced the short-hand notation 

r 2 (i- e )r(i + £ ) 
cr = r(i - 2e ) ' (35) 

Q 2 is an auxiliary variable which cancels in the full expression, /i is the scale 
introduced by dimensional regularization, ip^ denotes the colour-connected 
Born squared matrix elements, C(a;) is the colour charge of parton ai and 
the constant 7 (a;) gives the size of the virtual contributions to the diagonal 
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Altarelli-Parisi kernel P a ,/ai(0 

C(g) = N c ; 7 (.g) = UNc ~ m for gluons (36) 

o 

iV 2 - 1 3(N 2 - 1) 

C{q) = iiiVp 7(g) = AN for quarks , (37) 

finally V'ns^ represents the remaining finite terms- ._- 

The derivation of this result is rather simple ErEl First we observe that 
in axial gauge the collinear singularities come from the self energy corrections 
to the external lines □. For each helicity and colour sub-amplituds they are 
proportional to the Born term since the Altarelli-Parisi functions, P a / a (z), for 
diagonal splitting preserve helicity in the z — > 1 limit. Therefore, the collinear 
singularities of the one-loop amplitudes have the form 

4r=-(£) 2 E^ tree - ( 38 ) 

a 

There is a contribution for every external leg and the full contribution to -0( n ' 1 ) 
is easily obtained using eq. (|30|). 

The structure of multiple soft emission from hard processes in QED was 
investigated by Grammer and Yennieo. They have shown that the energetic 
electrons participating in a hard process receive an eikonal phase factor. In 
quantum chromodynamics, the situation is very similar except, that the eikonal 
factor is a-matrix equal to the path ordered product of the matrix- valued gluon 
fieldEJM 

For one soft gluon, the main result is very simple: it states that the singular 
contributions come from configurations where the soft gluons are attached to 
the external legs of the graphs. Therefore, the soft contribution can easily be 
calculated in terms of the Born amplitude. The insertion of a soft gluon that 
connects the external legs i and j has a twofold effect. First, after carrying 
out the loop integral and dropping singular terms corresponding to collinear 
configurations, we pick up the same eikonal factor as in QED f| 

9\\ 1 ( M 2 



^— it;*? 1-5- 1- (39) 



Secondly, the remaining part of the amplitude is the same as the Born am- 
plitude except that it gets rotated in the colour space by the insertion of the 

d Using unitarity this form of the soft factor can be confirmed by integrating the gluon 
momenta over the bremsstrahlung eikonal factor, as we shall see in the next section. 
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colour matrices appearing in the two vertices of the soft line connecting the 
hard lines i and j, therefore we have the replacement 



•^ Cl c 2 ..c„ ^ ^3 2-~i L c,c i CjCjC 1 ...c'...c'....c„ ' 



where Cj's denote colour indices for the external partons and t a , is the SU(3) 
generator matrix for the colour representation of line i, that is is (1/2)A° 7 - 
for an outgoing quark, — (1/2) A* ° for an outgoing antiquark, and ifaij for an 
outgoing gluon. For an incoming parton, the same formula can be used as 
long as we use the conjugate colour representations, — (1/2)A* J ° for a quark, 
(1/2) Xfj for an antiquark, and ifaij for a gluon. Finally, using the definition 
of V (n,1) we obtain 



ipij'°\{ai}i, n , {Pl}l,n) 



1 




PjUj{ai)u)(a2) 



2TZeal | >, >, >, ^^^^.^^^ |(41) 



This result of eq. (35) is scheme-dependent 0. In conventional dimensional 
regularization ^( 3 >°) and ipmn^ denote the Born and colour-connected Born 

(3 1) 

cross section in 4 — 2e dimensions and ip N 'g denotes the remaining final terms 
in this scheme. Finally, we note that to obtain the form of eq. (E3l one should 
expand the eikonal factor in e and apply the soft-colour identityT3 



3=1 

i¥=3 



4?'°\{Mi,ndPl}i,n) = 2C(a i )^ n < \{a l } hn ,{ Pl } ltn ). (42) 



6 Real contributions 

In this section I consider the limiting behaviour of the real contribution of the 
process (|24|) in the soft and collinear limit. I also give the local subtraction 
terms which render the real contributions integrable over the whole phase- 
space region. The discussion will be detailed on the soft contributions. In the 
case of the collinear limits I only summarise their most salient features. 
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6.1 Kinematics 



The four-momenta of the reaction ( p4| ) can be parameterised, for example, in 
terms of transverse momenta, rapidities and azimuthal angles 

Pi = "^(^1,0, 0,xi) , p% = ^-{x 2 ,0,0,-x 2 ) 

Pi = Px,i (ch yj, cos (fo, sin sh yi), ie{3,...,n}. (43) 

From energy and longitudinal momentum conservation we obtain for the mo- 
mentum fractions 

^ n 1 71 

xx = -p ^2pie Vi , x 2 = —^2 Pie -y* . (44) 

Considering singular limits we shall always assume the use of some suitable 
set of independent variables appropriate for the phase-space integration after 
imposition of four-momentum conservation. For the definition of measurable 
quantities, rapidity and transverse momentum variables are particularly con- 
venient since they are boost-invariant 1MB. In the evaluation of the soft and 
collinear limit it appears more convenient, however, to use energy-angle vari- 
ables 

= (Ei,sm<f)ism6i,cos(l>iSm.9i,cos9i) . (45) 



6.2 Soft subtraction terms and soft contributions 

The cross section of the real contribution (see eq. (pl|)) is constructed from 
products of three factors: function ip( n+1 '°\ the measurement function Sx.n-i 
and the phase-space factor d<j) n -i. Considering their soft limits let us assume 
that parton a k is soft (3 < k < n + 1). Its energy is denoted by E k and its 
angular correlations are controlled by the four- vector defined by the relation 

k^ = E k n^ k =E k (l,n k ). (46) 

The method that has been used in the previous section to calculate the soft 
limit of the virtual contribution also applies for the real emission, and 

we obtain 

47ra,c 



lim V (n+1 ' 0) ({ai}i,n+i;{Pi}i,n+i) = S ga k —^- Y] e lJ {p i ,p J1 n k ) x 

1>fr°\{aitfLil{pi&) (47) 
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where the list {uj}m m -i_„ denotes the same list as {vi} m _ m + n but Vk (m < k < 

m + n) is left out, ipij^ is the colour-correlated Born contribution defined by 
cq. ( |4l| ) in the previous section and e(pi,pj,nk) is the eikonal factor for real 
emissions: 

e(pi, Pj ,n k )= Pi ' Pj . (48) 

Pi ■ n k Pj ■ n k 

We note that e(pi,pj, nk) is independent of Ek but is dependent on the angular 
variables of the soft line. The colour-correlated Born terms ipij, however, are 
completely independent of the soft momenta. The soft limit of the measure- 
ment function is given by the requirement of infrared safety 

lim Sx, n -i({pih,n+i;X) = Sx,n~2({pi}t} n+ i;X) (49) 

Uih— >u 

and the phase-space factor behaves as 



lim i 

E k ^0 



i-l(Pl,P2 -> {Pl}3,n+l) = ——-d<f>l[k]d<j>n-2 (pi,P2 ~> {P^n+l) 



(50) 



where 



= ^ I T^rho = ^ I nrfC-SA ^ I d ^-2e (51) 



(2 7 r) d - 1 2p° ^ J 2(27r)( 3 - 2e ) 

is the phase-space integral over parton and the decomposition into energy 
and angular integrals is also shown^J. 

The local soft subtraction term for subtracting the soft singular behaviour 
in the Ek — * limit is given by the integrand of the expression below 

e [*] 

(soft, S ub) ( } = _K^Os_ y 

ai a 2 ,k Wl'«> > n _ I 2lT 



E 

r 1 [*] 



. 1 



87r e(pi,pj,n k ) -^2®{Ek > E c )d4>i[k] 



El 



^ - 2 y d<Pn~2{pi,P2 ~> to>3 fc L+l) (52) 



e We note that in the soft limit the condition imposed by the delta function of momentum 
conservation is different from the original one; therefore our notation is not completely 
precise. In the soft limit some of the components of the momenta {pi}^^ i j will be different 

from the original ones and instead {p;}^,).! we should write lim^;,.— >o{P!}3 C | l _i_i • We tacitly 
assume that in the soft limit pi denotes its limiting value. It is uniquely defined after choosing 
the independent set of variables in the phase space-integral drf> n —i. 



1G 



In eqs. ( j3l| ) and (p2|), wc can use the same independent integration variables, 
therefore, by adding eq. ( |52| ) to the integrand of ([H]) we subtract its singular 
regions defined by the E/~ — > limit. Again, to ensure that the subtraction 
does not change the value of the original expression what is subtracted has to 
be added back on but in an integrated form where the singular terms are cal- 
culated analytically. This is achieved by carrying out the phase-space integral 
d(j)i [k] . The expression in the square bracket in eq. ([32]) defines the soft integral 



X^ ott (z„-, M /£; c ,e)-87rV 



Pk 



Pi ■ Pi 



2p° k (2ir) d 1 pi -pkPjPk 



-Q{E k > E c 



(53) 



where E c is a cut-off value on the energy of the particle which is allowed to be 
soft. We also indicated the remaining dependences of the integral; the angular 
variable z a b is defined as 



z ab = cos V ab 



PaPb 



Wa\ 



Definition 


Answer 


T / \ Oi_2e fl j n f 71 " it (sin 9 sin <b)~ 2e 

J ( e ) = nl.ll I-i dcosd Io W i-cosfl 

_ rl (1-z 2 )- 

J-l 1-2 


-i + 21n2 
+e(Li 2 (l)-21n 2 2) 




7T 


If(Zij , Ztk) — J Q d<t>k i~ ZtjZik - s i n e tj sine ifc cos0 fc 


\Zik~Zij | 


jW(z ij ) = ±f 1 _ 1 dz ij f*d ( l> k lx 




J W (zij) = \ J-i dzn ft d(j> k In sin 2 <j> 1% 


-In2(l + Zy)hii^ 


J (z \zij) = i/^^ lnsin 2 ^/;^ l* 


i In 2 2(1 - Zij) - 2 In 2 2 
+Li 2 (l) - Li 2 (^p) 



Table 2: List of angular integrals 

It is convenient to decompose the eikonal factor into a term which has one 
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collincar singularity and terms which arc finite 



where 



(Pi -Pk)(Pj -Pk) 



E\ 



2 [K Z iji Z ik, z ]k) + {i ^ j) 



1 - Zik 



H - I {z%j i Zik-> Zjk) 



and 



l 



(1 - Z ifc )(l - Zjfc) 



1 



Zik 



1 - ^ 7 fe 



(54) 
(55) 

(56) 



The integral over the energy is trivial and the angular integrals have to be 
calculated only up to order e so one obtains 



j soft 
1 ij 



*<^>-(-i) 

J{e) + J (0) (%■)(! " e In 4) - eJ (z) (%) + J^H^)] (57) 



where cr was defined in eq. @. The integrals J^, J<» and jW are listed in 
Table (||). Inserting their values into eq.(^7|) and expanding in e we get both 
the singular and finite parts 



j soft 

a 

j soft, sing 
ij 



j soft, sing 



r soft. fir. 



1 1 ]u ZPiPi 
e 2 e Q 2 



1 EiEj 
- In — tt- 

e E 2 



(58) 
(59) 



The explicit farm of the finite part is not interesting for us here; it can be 
found in ref. Ell. We can use completely covariant notation by replacing the 
energies with the covariant expressions 



Ei = (Vpi)/yfi , where V = p\ + p% = (y^, 0, 0, 0) 



(60) 



where pa and ps are the four momenta of the incoming hadrons. 

If the sum over flavour is carried out the result becomes independent of 
the label of the soft line. Ttherefore every leg in the final state gives the same 
contributions and the factor l/(n — 1) in eq. ( 521 ) gets cancelled. As a result, 
the soft contribution of the real leading order process with n + 1 partons can be 
written in the form of the virtual contribution of the corresponding n-parton 
process. It can be obtained form the right hand side of eq. (p2) by multiplying 
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it with n — 1, changing its sign and inserting the integrated value of Jf° 
(eq.|59|): 

dai n + 2 1 ^(p 1 ,p 2 ;X) = ]T ^ n+hsoit) S x ^ 2 ({ Pl } 3 , n ;X)x 

{ai} 3 , n 

d<f> n -2(pi,P2;{pi}a,n) (61) 

where 

{t^ + ^§-} ^ ({a l}hn , { Pl } hn ) 

+ ^S +1 ' SOh \Wlh,n,{PlKn) (62) 

The first term is the soft-collinear singularity and it cancels the soft-singular 
terms appearing in the virtual corrections. The third term is the soft contribu- 
tion proportional to the colour correlated Born-term and again it cancels the 
corresponding terms in the virtual contributions. The second term comes from 
the collinear singularities of the eikonal factors and is cancelled by the direct 
singular collinear contributions. The last term is finite and can be evaluated 
in four dimensions. 

6.3 Collinear subtraction terms and collinear contributions 

In the previous section we constructed the local subtraction term for the soft 
singular region. We also demonstrated that the soft singularities of the virtual 
corrections are cancelled by the soft contributions of the real corrections after 
the integrals over the energy and angular variables of the soft line are carried 
out. 

The subtraction and addition procedure can also be applied to the singular 
collinear regions. Similarly to the soft case, one finds again simple limiting be- 
haviours for the ^-functions, the measurement functions and the phase space. 
We shall discuss only the collinear limit of the %p function. The reader can find 
further details in refs.li34§y. 
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From the simple behaviour of the helicity amplitudes in the collinear limit 
we obtain for the ip function the limiting behaviour 



K ' Pj 
4iras 



+ -Qa t S(a^)*(z)i> in ' Q) (Pl:2; .., Pi , .., Pj , ..) , (63) 

Pi-P 3 

where z denotes the momentum fraction defined through the equations 

Pi = zp P , kj = (1 - z)p P , (64) 

Pab(z) denotes the standard Altarelli-Parisi splitting functions and Q a b*(z)'s 
are some new universal functions which control the azimuthal angle behaviour 
of the collinear limit. Assuming that the collinear particles belong to the final 
state 



-4C A z(l-z), (65) 

4T F z(l-z), (66) 

0, (67) 

0. (68) 

In these equations, the * symbol over the flavour of the particle that eventually 
splits reminds that this particle is off-shell. In principle, this notation should 
be extended also to the Altarelli-Parisi splitting kernels, but at the leading 
order P^* = P a *b, and therefore there is no need to keep track of the off- 
shell particle; the ip function is constructed from the helicity amplitudes of 
n-parton processes but with some linear dependence on the azimuthal angle of 
the quasi-collinear configuration. The ip functions therefore are as simple as 
the Born-terms. They are important in constructing a correct local subtraction 
term. However, upon integrating over the azimuthal angle of the collinear 
momenta their contributions vanish and they do not appear in the integrated 
collinear contributions. It is an interesting simplifying feature of the five parton 
processes that the ip functions vanish identically 

Using crossing properties of the splitting functions and of the ip functions 
similar relations remain valid also for initial collinear singularities. 

The remaining procedure, in principle, is the same what we used in the 
soft case although it is somewhat more tedious. The properties of the S- 
functions and the phase-space again ensure that we can easily construct the 



Qgg*( z ) = 

Qqg* ( z ) = 

Qgq* (z) = 

Qqq*{z) = 
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local collinear counter terms and in the collinear contributions we can obtain 
the singular terms analytically. For further details the reader should consult 
the original literatureEZrE£l. 

Exercise: Calculate the inclusive one-jet cross section da/dEjdcosQj for 
e + e~ annihilation following the general algorithm explained above. 



7 Conclusion 



In this lecture I described the methods and techniques to get differential cross- 
section formulae from NLO singular scattering amplitudes. These are free from 
singular terms well defined in all integration regions and therefore suitable for 
numerical evaluation, The ingredients of the methods are the collinear counter 
terms, the local subtractions terms, the virtual, the soft and the collinear 
contributions. Their sum defines the finite hard scattering cross-section 



j-hard j_born , j_virt , 7 soft , j_(coll, initial) , j_(coll. final) 

"Oajas = d V ai a 2 + rfo a 1 a 2 +« cr a 1 a 2 + dcr a 1 a 2 + rfcr a 1 a 2 

I ^^.(coll, counter) _|_ ^^.(real, subtracted) (69) 



where the last term has the kinematics of an n + 1 parton process while all 
the other terms have the kinematics of an n parton process. In addition, 
the evaluation of some physical quantity requires the explicit construction of 
the corresponding measurement functions. We "P^ that in the case of jet- 
production this is a non- trivial exercise (see refsEZJE3rEj). 

We have seen that although the methods are conceptually simple, their 



implementation is non-trivial. , ^, pr , PT 
comprehensive documentation 
reading. 



several papers attempted to give a 
which are recommended for further 
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